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Abstract 

PHANTOM is a tree level Monte Carlo for six parton final states at proton-proton, 
proton-antiproton and electron-positron collider at 0{a%j^) and C(Q;^J^^a|) includ- 
ing possible interferences between the two sets of diagrams. This comprehends all 
purely electroweak contributions as well as all contributions with one virtual or two 
external gluons. It can generate unweighted events for any set of processes and it is 
interfaced to parton shower and hadronization packages via the latest Les Houches 
Accord protocol. It can be used to analyze the physics of boson boson scattering, 
Higgs boson production in boson boson fusion, tt and three boson production. 
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PACS: 12.15.-y, 11.15.Ex 

Classification: 11.1 General, High Energy Physics and Computing 

External routines /libraries: LHAPDF (Les Houches Accord PDF Interface), CIRCE 

(beamstrahlung for e"'"e~ ILC collider). 

Nature of problem: 

Six fermion final state processes have become important with the increase of col- 
lider energies and are essential for the study of top, higgs and electroweak symmetry 
breaking physics at high energy colliders. Since thousands of Feynman diagrams con- 
tribute in a single process and events corresponding to hundreds of different final 
states need to be generated, a fast and stable calculation is needed. 
Solution method: 

PHANTOM is a tree level Monte Carlo for six parton final states at proton-proton, 
proton-antiproton and electron-positron collider. It computes all amplitudes at 
^{'^%m) ©(a^ji^al) including possible interferences between the two sets of di- 
agrams. The matrix elements are computed with the helicity formalism implemented 
in the program PHACT [1]. The integration makes use of an iterative-adaptive mul- 
tichannel method which, relying on adaptivity, allows to use only few channels per 
process. Unweighted event generation can be performed for any set of processes 
and it is interfaced to parton shower and hadronization packages via the latest Les 
Houches Accord protocol. 
Restrictions: 

All Feynman diagrams are computed al LO. 

Unusual features: 

PHANTOM is written in Fortran 77 but it makes uses of structures. The g77 com- 
piler cannot compile it as it does not recognize the structures. The vintel, Portland 
Group, True64 HP fortranTT or fortranQO compilers have been tested and can be 
used. 

Typical running time: 

A few hours for a cross section integration of one process at per mille accuracy. One 

hour for one thousand unweighted events. 
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ph/9403244]. A. Ballestrero, PHACT 1.0- Program for Helicity Amplitudes Cal- 
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LONG WRITE-UP 



1 Introduction 

Monte Carlo event generators are essential tools for the comparison of theory 
and experiment in High Energy Physics. In this paper we present PHANTOM a 
new event generator which is capable of simulating any set of reactions with 
six partons in the final state at pp, pp and e+e" colliders at 0(a|j^^) and 
0{a'^j^al) including possible interferences between the two sets of diagrams. 
This includes all purely electroweak contributions as well as all contributions 
with one virtual or two external gluons. The relevance of six fcrmion final states 
both at the LHC and the ILC is well known and will be further discussed in 
Sect. 2. The signal of processes like Higgs boson production in vector boson 
fusion with a Higgs decaying in two bosons, boson boson scattering and 3- 
vector boson production are all described by 6 fermion final states at 0{a^^). 
Top-antitop pair production is also described by a six fermion final state and its 
signal is at 0{a^j^) for the ILC but at hadron colliders the main contribution 
comes from 0{a'^^a1). 

As far as the irreducible backgrounds to the above mentioned signal is con- 
cerned, 0{aj^j^ja1) reactions encompass the full lowest order QCD background 
for six-parton studies of final states with four leptons at hadron colliders. An 
example is VV fusion with fully leptonic decays of the intermediate vector 
bosons. On the contrary, semileptonic channels get additional contributions 
from 0{a%^a^g) processes which are responsible for the V+Ajets background. 
The latter must be included as well in a complete analysis and can be covered 
by other Monte Carlo event generators, such as AlpGen [1] or MadEvent [2]. 
It should be noticed nevertheless that these contributions have quite different 
kinematical features with respect to the scattering signature, therefore they 
are expected to be suppressed by means of appropriate selection cuts. 

At e+e^ colliders there is no QCD background for six and four lepton final 
states and 0{a^^al) reactions describe completely the QCD background to 
final states with two leptons. Only for fully hadronic final states the set of 
reactions presently available in PHANTOM needs to be complemented by other 
generators. 

PHANTOM employs exact tree matrix elements. There are issues which cannot be 
tackled without a complete calculation, like exact spin correlations between the 
decays of different heavy particles, the effect of the non resonant background, 
the relevance of the offshcllness of boson decays, the question of interferences 
between different subamplitudes. Without a full six parton computation it will 
be impossible to determine the accuracy of approximate estimates. In Ref. [3] 
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the complete calculation for PP ^jfJ-^fj. at 0{a%f^) has been compared at 
length with a production times decay approach, showing differences of the 
order of 10-20% in some important regions of phase space. The reliability 
of the Equivalent Vector Boson Approximation (EVBA) [4] in the context of 
vector boson scattering has been critically examined in [5]. 

PHANTOM is an example of a dedicated event generator which describes a pre- 
defined set of reactions striving for maximum speed and efficiency. Other re- 
cent examples of dedicated programs for LHC physics are Toprex [6], which 
provides the matrix elements for several reactions related to top production, 
Alpgen [1] and GrOppa [7]. Alpgen can simulate a large number of processes 
with electroweak bosons and heavy quarks plus jets. When more than one 
electroweak boson are present in the final state they are considered to be on 
shell. GrSppa provides event generators for V (W or Z)+ jets (< 4 jets), VV 
+ jets (< 2 jets) and QCD multi-jet (< 4 jets) processes. 

Dedicated generators aimed at future e^e~ colliders are Lusifer [8], which 
covers all reactions with six massless fermions in the final state at O^a^j^^^) 
and 0{ai^al), SIXFAP [9], which covers all reactions with six fermions in the 
final state at 0{a%j^) including mass effects, SIXPHACT [10], which was among 
the first MC to compute six fermion final states but was limited to processes 
with one final neutrino, eett6f [11], which simulates only processes related 
to top-antitop production and SIXRAD[12], which deals with all six jet final 
states at 0{a1,^a\). 

The Monte Carlo's just mentioned show that six fermion physics has been 
already investigated at the ILC since a long time and they have been used 
for several phenomenological studies. PHANTOM is the first complete dedicated 
Monte Carlo which extends this kind of physics to hadron colliders. 

The complementary approach is given by multi-purpose programs for the au- 
tomatic generation of any user-specified parton level process. The following 
codes for multi-parton production are available: Amegic-Sherpa [13], CompHEP 
[14], Grace [15], Madevent [2], Phegas & Helac [16], O'Mega & Whizard [17]. 
CompHEP is limited to processes with at most eight external particles but since 
it computes matrix elements squared instead of helicity matrix elements it 
is extremely slow when a large number of particles is involved. Grace is a 
framework for generating single process matrix elements which can be in- 
terfaced to integration and event generator modules. The remaining pack- 
ages, Amegic-Sherpa, Madevent, Phegas & Helac and O'Mega & Whizard 
are complete event generators which automatically generate the amplitudes, 
produce the mappings for integration over phase space, compute cross sections 
and generate unweighted events. While processes with six particles in the fi- 
nal states can be tackled with these generators, see for instance [18], they are 
at the limit of their practical capabilities when the unweighted generation of 
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several hundreds of this kind of processes is involved, as it is needed for LHC 
physics. 

PHANTOM profits from the experience obtained with PHASE [19] which could 
only simulate all electroweak processes of the type PP — > 4g + lui at the 
LHC. The two codes share a number of key features. The matrix elements 
are computed with the use of the modular helicity formalism of Refs. [20,21] 
which is well suited to compute in a fast and compact way parts of diagrams 
of increasing size, and recombine them later to obtain the final set. The in- 
tegration strategy merges the multichannel approach [22] with the adaptivity 
of VEGAS [23]. This results in generators which adapt to different kinemati- 
cal cuts and peaks with good efficiency, using only few channels per process. 
Both codes employ the one-shot method developed for WPHACT [24], and used 
for four-fermion data analyses at LEP2. In this running mode, all processes 
are simultaneously generated in the correct relative proportion for any set of 
experimental cuts. 

The paper is organized as follows: in Sect. 2 we examine the physical content of 
the processes which can be described by PHANTOM. Then Sect. 3 describes the 
features of the program, the way amplitudes are calculated, phase spaces are 
implemented and integration is performed. The following section is dedicated 
to the modes of operating the program and finally Sect. 5 reports an example 
of the phenomenological studies which can be performed using PHANTOM. 



2 Physical processes 

Processes with six partons in the final state will be central to the physics 
program at next generation accelerators, the LHC [25,26] and the ILC [27]. 
They include Higgs boson production in vector boson fusion followed by Higgs 
decay to WW and ZZ, boson boson scattering processes, top-antitop pair 
production, three vector boson production. 

The Standard Model (SM) provides the simplest and most economical ex- 
planation of Electro- Weak Symmetry Breaking (EWSB) in terms of a single 
Higgs doublet. The search for the Higgs sector and its investigation will have 
the highest priority for all LHC experiments [28-30]. Detailed reviews and 
extensive bibliographies can be found in Refs. [31-33]. In the SM the Higgs is 
essential to the renormalizability of the theory and is also crucial to ensure 
that perturbative unitarity bounds are not violated in high energy reactions. 
In particular, if no Higgs is present in the spectrum, longitudinally polarized 
vector bosons interact strongly at high energy, violating perturbative unitarity 
at about one TeV [34]. If, on the contrary, a relatively light Higgs exists then 
they are weakly coupled at all energies. These processes have been scrutinized 
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since a long time, going from the pioneering works in [35,36], which address 
boson boson scattering on a general ground, to the more recent papers in 
[37,38] focused on the extraction of signals of vector boson scattering at the 
LHC. In the last few years QCD corrections to boson boson production via 
vector boson fusion [39] and to three vector boson production [40] at the LHC 
have been computed. The size of the corrections depends quite strongly on 
the particular process under consideration: while in the boson fusion case the 
corrections are below 10%, for the three boson production case they are in the 
70 to 100% range. 

Top pair production with its large cross section and large signal to background 
ratio will play a central role since the early days of data taking at the LHC. 
The measurement accuracy will quickly be dominated by systematic effects, 
in particular by 6-jet energy scale uncertainties, while statistical errors will 
be neghgible. The wealth of available channels will provide feedback on the 
detector performance and assist in detector calibration [25,26]. A total error 
on the top mass of about 2 CeV is expected to be feasible. Exploiting the 
high statistics it will be possible to search for new physics signatures in top 
production and decay. Spin correlations will be measured to a precision of a 
few percent. Finally, a good understanding of top physics will be essential since 
top production will be a major source of background for all investigations of 
processes involving high pt H^'s, like boson boson scattering and many new 
physics searches. 

PHANTOM is intended in particular for studies of Higgs, boson boson scattering 
and top physics in six fermion final states. It can compute all processes of this 
kind and does not make use of any production times decay (narrow width) or 
equivalent vector boson approximation. This implies that every single diagram 
contributing to a definite final state is included and not only the resonant 
ones. So for instance if one studies a vector boson fusion process not only the 
diagrams reported in Fig.l will be considered but also, when appropriate, all 
other diagrams of the types described in Figs. 2,3,4. In the same way when 
studying top antitop production not only the first two diagrams on the left of 
Fig. 4 will be considered but all the others schematically described in the rest 
of Fig. 4 and also in the other 3 figures. Moreover all diagrams which do not 
have two resonant bosons in the final state contribute as well. We have already 
discussed in some particular cases [10,3,5] the differences between a complete 
calculation and an approximated one. These depend on the particular physical 
study at hand and on the applied cuts. But a complete calculation is in any 
case an important tool to quantify the possible discrepancies. 

Let us explain with one example which are the processes to be included in a 
practical case. Suppose one wants to study a Higgs produced in boson boson 
fusion decaying to two W's with one of the W^s decaying leptonically. In such 
a case one will experimentally search for a signal of four jets plus a lepton and 
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Fig. 1. Vector boson fusion processes. 




Fig. 2. Non fusion and non doubly resonant two vector boson production. Three 
vector boson production. 




Fig. 3. Higgs boson production via vector boson fusion and Higgsstrahlung. 




Fig. 4. Examples of contributions to the QCD irreducible background: tt production 
(a,b) and VV + 2 jets (c,d) 

missing transverse momentum. The corresponding physical processes that will 
have to be generated with PHANTOM are all those with a lepton and a neutrino 
and 4 quarks in the final state at + (9(Q;|j,^a|) as well as those with 

one or two outgoing gluons and three or two quarks in addition to the leptons. 
The number of these processes amounts to a few hundreds. The unweighted 
generation of all of them will be performed simultaneously with PHANTOM. 

Analogous considerations hold of course for top antitop or boson boson scat- 
tering studies. One cannot separate the different physical processes just con- 
sidering the final states, and the strategy is that of a full calculation. With 
a complete sample one will thereafter evidentiate the physics of interest with 
appropriate cuts. 
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3 Program Description 



We give in this section a description of the main features of the PHANTOM Monte 
Carlo. The routine hst and an example of input file for running the program 
are presented in appendices A and B respectively. A User Guide can be found 
in the PHANTOM distribution, with more details on how to run the program and 
how to prepare the input files. 



3. 1 General features 

As already said, PHANTOM can generate unweighted event samples for any set 
of reactions with six partons in the final state at pp, pp and e'^e~ colliders 
at 0{a%j^) and 0{a\^al). Possible interferences between the two sets of dia- 
grams can be taken into account. Alternatively, the user can simulate purely 
electroweak and 0{a%^^a'j?) processes separately. All finite width effects are 
taken into account as well as all correlations between the decay of any un- 
stable particle which appear in intermediate states. In most cases, with the 
notable exception of tt production, the purely electroweak contribution con- 
tains the physical processes we wish to investigate, while the contributions 
with one virtual or two external gluons represent a background. In particular 
0{a%j^a1) processes completely describe Wjj productions at hadron colliders 
and V4j' at e'^e~ machines. 

The program proceeds in two steps which are performed separately. In the first 
phase each of the reactions which need to appear in the final event sample 
is integrated over and the corresponding cross section is calculated. At the 
same time the density of points in the unit cube employed for the sampling is 
recursively modified in order to match the actual behaviour of the integrand 
and the integrand maximum is searched for. The information gathered in this 
step is stored in the form of an integration grid. In the final stage all the 
grids and maxima are used to produce the unweighted event sample. The 
separation between the two phases allows for additional flexibility: the grids 
produced with a specified set of acceptance cuts can be used to produce event 
samples with any set of more restrictive selection cuts or for a smaller number 
of reactions without repeating the integration step. Additional processes can 
easily be introduced. 

The number of reactions which contribute to a given final state, say for in- 
stance 4j7z/, can be very large, particularly at hadron colliders. In our approach 
it is crucial to keep the number of reactions for which the individual cross sec- 
tion needs to be computed to a minimum. By making use of symmetries, the 
task of computing all relevant processes can be substantially simplified. Indeed 
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all reactions which differ by charge conjugation and by the symmetry between 
the first and second quark families can be described by the same matrix ele- 
ment, provided one take a diagonal CKM matrix and ignores the light quark 
masses. This extends to processes which are related by a permutation of the 
lepton families. At most they differ by the PDF's convolution and possibly by 
a parity transformation, which can be accounted for automatically. In the in- 
tegration step the user can sum over all processes which are related by charge 
and family symmetries setting the flag i_ccf am. 



3.2 Amplitudes 

The matrix elements are computed with the helicity formalism of Refs.[20] 
and the semi-automatic method described in Ref.[21]. In our experience the 
resulting amplitudes are much faster than those obtained by completely auto- 
matic procedures. Even though wc have not made a detailed comparison with 
all programs and for all kind of amplitudes, we happened for instance to find 
differences in CPU time of one order of magnitude or more with Madgraph, in 
agreement with the statements in Ref . [39] . 

In our calculation wc first of all take advantage of the well known fact that all 
processes which share the same total particle content, with all partons taken 
to be outgoing, can be described by a single master amplitude. 

The calculation can be further simplified examining more closely the full set 
of Feynman diagrams. Since we assume a diagonal CKM matrix, the fcrmion- 
antifcrmion pair appearing at the end of a given fcrmion line cither are of the 
same flavour, if only neutral vector insertion or an even number of W inser- 
tions take place along the line, or reconstruct an isospin doublet if an odd 
number of W insertions takes place. Therefore, limiting ourselves to reactions 
in which all external particles are fermions it becomes quickly apparent that 
in some processes, fermions can be paired only into neutral couples (we in- 
dicate them as 4Z), e.g. uc — > uciijiee, while in others only charged couples 
can be formed (4W), e.g. ud — > cs^VuVeB. In other cases they can form two 
charged and two neutral couples (2Z2W), e.g. uu hb^V^Vee. In all other 
cases the full set of diagrams splits into a sum of the three basic gauge invari- 
ant sets mentioned above, and can be described as 2Z2W+4Z,4W+2Z2W or 
4W+2Z2W-I-4Z, e.g. ud — > udeV^u^e. Similarly, all reactions with two external 
gluons and six fermions can be computed from two basic amplitudes which, 
with obvious notation, we call 2g2WZ and 2g3Z. Prom the above discussion 
it becomes clear why in PHANTOM we have eight master amplitudes: three pure 
electroweak with eight external fermions (4Z, 4W, 2Z2W) at 0{a%j^j)^ the same 
three at + O^a^^a^g) and the two with two external gluons (2g2WZ 

and 2g3Z). We take moreover advantage of the fact that the diagrams in which 



9 



two or more of the external particles are identical can be split in sets which 
just differ for the exchange of the identical particles. Therefore the master 
amplitude corresponds to the basic diagrams that one would have if all ex- 
ternal particles were different. In PHANTOM the full amplitude of the process 
at hand is then computed by the routines ampem.f , ampSfqcd.f , amp2g.f 
which call all the necessary master amplitudes with the appropriate order of 
the momenta and call coleval.f or colevalew.f to evaluate the color fac- 
tors. The amplitude routines compute separately the sum of the diagrammatic 
contributions for every possible color configuration. The probability used for 
choosing the color fiow of each event, as specified by the Les Houches Ac- 
cord File Format, is taken to be proportional to the modulus square of the 
amplitude of each conguration. 

The range of reactions which are available in PHANTOM is substantially larger 
than in PHASE. All 0{a^j^al) contributions with one virtual or two external 

gluons and the 4Z amplitude were not previously available, only two master 
amplitudes were needed and general routines to compose the different ampli- 
tudes and evaluate the color factors were not necessary. 

3.3 Master amplitudes computation 

In this section we want to show briefly how the computation of master ampli- 
tudes is organized. We will flrst consider a master amplitude for eight external 
fermions at with one gluon exchange and then the one for 2g2WZ. 

We refer to Ref.[19] for a description of the kind of strategy we use for 0{a^j^) 
eight fermion master amplitudes. 

In the flrst case the possible ways in which the gluon can be exchanged are 
depicted in Fig. 5 




Fig. 5. Gluon exchange diagrams 

This immediately shows that it is convenient to perform the calculation via the 
subprocesses of Fig. 6. Starting from the left, the first two correspond to the 
decay of a virtual boson or gluon to four fermions, the third and the fourth are 
building blocks in which two external fermions annihilate in a virtual boson 
or gluon which is inserted in the middle of a fermion line. At least one of the 
two ends of the fermion line does not correspond to an external particle. The 
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other parts of the fermion hne and its attachments can be easily combined 
with such building blocks in the method we use [20,21]. 



Fig. 6. Gluon exchange subdiagrams 

For the master amplitude of 2g2WZ we notice first of all that the possible 
electroweak configurations are those of Fig. 7. To obtain all possible diagrams 
of the master amplitude it is sufficient to attach in all possible ways the two 
external gluons. 



Z 

W 

w 



z 

w 
w 



x2 



Tzih 

x2 



Z 

w 
w 




Fig. 7. Electroweak configurations of 2g2WZ 

This leads to the twelve color structures of Fig. 8 where the first two starting 
from left correspond to both gluons attached to the same quark line in all 
possible ways and the second two to the two gluons attached to two different 
quark lines. It is superfluous to remind that the configuration with the two 
external gluons forming a triple vertex with a gluon propagator attached to 
the quark line can always be reduced to a linear combination of configurations 
a) and b). 
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Fig. 8. Color structures of 2g2WZ 

Also for 2g2WZ we proceed evaluating subprocesses in order to avoid com- 
puting twice a subdiagram appearing in several different feynman diagrams. 

Therefore one starts with subdiagrams of increasing complexity: virtual boson 
decaying to two fermions, to two fermions and one gluon and to two fermions 
and two gluons. Then those corresponding to one virtual fermion decaying to 
three fermions or three fermions + 1 gluon, etc. Putting together all these 
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pieces one finally arrives at computing all groups of Feynman diagrams be- 
longing to the the same color structure, which is precisely the purpose of the 
routines computing the master amplitudes. 

3.4 Parameters and Parton Distribution Functions 

Standard model parameters are defined in the routine coupling. f. In our 
notation, rmw, rmz, rmt, and rmb are the W, Z, top and bottom masses re- 
spectively. The total W and Z widths are given by gamw and gamz. Higgs and 
top widths are computed in the same routine by standard formulas. 

PHANTOM employs the G^u-scheme defined by the input set: Mw, Mz and Gp- 
According to this scheme, the calculated parameters are 



sin^ 1 - {M^/Mif aem{Mw) = — GpM^ sin^ 9^ 

TT 

where 9^ is the weak mixing angle, and the electromagnetic fine structure 
constant. 

Wc introduce the decay width in the propagators of unstable particles using 
the fixed-width scheme. For the vector boson propagators, in the unitary gauge 
we work in, it consists in replacing with M^—iMF both in the denominator 
and in the p^p^ term. This scheme preserves U(l) gauge invariance at the 
price of introducing unphysical widths for space-like vector bosons. A general 
discussion of the fixed-width scheme as well as of alternative approaches to 
introducing decay widths in scattering amplitudes can be found in the papers 
in Ref. [43]. 

For the PDF's we use the Les Houchcs Accord distribution [44], which can be 
downloaded from http://projects.hepforge.org/lhapdf/ . 

3.5 Iterative-Adaptive multichannel 

For completeness we describe in this section the iterative and adaptive multi- 
channel technique which PHANTOM employs for integration. In our opinion the 
ability to adapt is the overriding consideration for multidimensional integrals 
of discontinuous and sharply peaked functions. An individual process can con- 
tain hundreds of diagrams and several unstable massive particles can appear at 
intermediate stages. The resonant peaking structure of the amplitude is there- 
fore generally very rich. As a consequence the 16-dimensional phase space has 
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multiple combinations of non-trivial kinematical regions corresponding to en- 
hancements of the matrix element which need to be smoothed out in order to 
achieve a good convergence of the integration process, which in turn results 
in good efficiency of event generation. 

We have merged the multichannel method [22] and the adaptive approach a la 
VEGAS [23]. An algorithm based on a similar philosophy has been proposed in 
[41]. In the multichannel approach, mappings into phase-space variables are 
chosen in such a way that the corresponding Jacobians cancel the peaks of 
the differential cross section. These mappings are not in general unique. One 
normally needs several different phase-space parametrizations, called channels, 
one for each possible peaking structure of the amplitude. This gives rise to a 
huge combinatorics which requires a correspondingly large number of channels. 
Number and type of these mappings must be fixed a priori, before starting the 
integration. The multichannel method thus requires a guess on the behavior 
of the integrand function. It indeed relies on the expectation that the selected 
set of channels, properly weighted [42], is able to describe reasonably well the 
amplitude. As no adaptivity is provided (apart from the freedom to vary the 
relative weight of the different channels), neglecting even one channel might 
worsen considerably the convergence of the integral. 

The criteria of adaptive integration as performed by VEGAS are rather differ- 
ent. This approach bases its strenght on the ability to deal automatically with 
totally unknown integrands. By employing an iterative method, it acquires 
knowledge about the integrand during integration, and adapts consequently 
its phase-space grid in order to concentrate the function evaluations in those 
regions where it peaks more. VEGAS divides the N-dimensional space in hyper- 
cubes, and scans the integrand along the axes. For a good convergence of the 
integral, it thus requires amplitude peaks to be aligned with the axes them- 
selves. The problem can be easily solved if one set of phase-space variables 
is sufficient to describe the full amplitude peaking structure. In this case, the 
alignement can always be obtained by an appropriate variable transformation. 
The method becomes inefficient when it is impossible to align all enhancements 
with a single transformation. 

In order to reduce the number of separate channels one has to consider in 
the multichannel, we use multi-mapping. In general multiple peaks can appear 
in the same variable, together with long non-resonant tails which extend far 
away from the peaks. This latter case is more and more severe as the collider 
energy increases. For instance the mass of a neutral fermion pair // would 
typically resonate at both the Z and Higgs mass. The two peaks and the three 
non resonant regions can be separately mapped into a five zone partition of 
the basic interval [0,1] of the integration variable. While multi-mapping is 
extremely useful to improve the convergence of VEGAS integration within a 
single phase-space parametrization, in general several such parametrizations 
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are needed. In this case, one has to introduce different channels (in standard 
multichannel language) with their proper multi-mapping. Each channel defines 
a non-uniform probability density which describes a specific class of 

amplitude peaks. One can then write 



N „ N 

Y^ocij dxi f'{G;\xi)) =Y.aJi (1) 



i=l " 1=1 



ttj being the so called weight of the z-th channel (,Xj = G'i($)), and /'($) the 
smoothed integrand. The aj quantify the relevance of the different peaking 
structures of the amphtude. Owing to the very poor knowledge of the inte- 
grand, it is rather difficult to guess these values a priori. Usually, they are 
computed and optimized during the integration run. In the iterative- adaptive 
multichannel^ the integral in Eq.(l) splits in distinct contributions. The 
presence of identical final-state particles increases the possible list of reso- 
nant structures. In order to keep the number of separate integration runs 
manageable, we include all jacobians generated by particle exchange in the 
denominator of Eq.(l), while exploiting the freedom to relabel the momenta 
to regroup all integration runs related by particle exchange to a single one. 



The integration proceeds through two steps: the first one which is called ther- 
malization collects preliminary information about the integrand and fixes the 
weight of the different channels. In every thermalizing iteration, all channels 
are independently integrated for some set of Oj. At the end of each iteration, 
a new set of phase-space grids (one for each channel), and an improved set of 
ttj are computed. The criteria for weight optimization we adopt is 



I 



(2) 



where Jj is the i-channel integral. The new set of and grids are then used in 
the next iteration. The procedure is repeated until a good stability of the a, is 
reached. Weights smaller than 1 x 10~^ after the third thermalization iteration 
are set to zero and the corresponding channels are discarded. If any number of 
channels is dropped two additional thermalization iterations are performed. In 
the second step no further change of the weights is allowed, the total integral 
for each channel is evaluated an a corresponding phase space grid is produced. 
These grids are then stored and made available for event generation. 
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3.6 Phase space 



While the number of possible combinations of resonances is very large, from a 
topological point of view the structure is much simpler. Apart from the weak 
vector bosons the only unstable particles in the SM which are relevant for LHC 
or ILC physics arc the Higgs boson and the top quark. While the weak bosons 
at tree level have only two body decays, the top quark decays to hW and then 
to a three body system. For the Higgs boson the situation is more complex and 
the ratio of the various branching ratios and therefore their phenomenological 
relevance depends crucially on the Higgs mass. The Higgs decay channels 
however fall into two groups: two body decays like 77, hb or r+r^ or four 
body decays through W^W^ or ZZ. It is therefore possible to construct a 
limited number of phase space parametrizations each of which, varying the 
masses and widths of the particles which appear in the intermediate states, 
can describe several channels in the multichannel sense. 

In addition to s-channel resonances, t-channel enhancement can also be present 
At the LHC they are regulated by minimum pt requirements. At the ILC the 
issue is much more relevant because processes with electrons lost in the beam 
pipe can be an important background to charged current reactions. 

Since we do not keep track of individual diagrams, the possible resonant struc- 
tures must be deduced from the external particle content only. This is per- 
formed in the proc.f routine. Notice that the resonant structure depends on 
the pcrturbative order: qq —>■ qqWW can have a if — > WW resonance at 
0{a%^) but not at 0{ai^al). 

The criteria we adopt to automatically define number and type of channels 
needed for a given process are the following: for each process we consider all 
particles as outgoing and then examine all sets of four fifj pairs which can be 
constructed, retaining only those in which all pairs can be identified with the 
decay products oi aW or Z boson. This produces the complete set of possible 
two-particle enhancements. As a second step, for each selected set of pairs we 
search for triplets which reconstruct bW'^ and bW~ states and W^W~ and ZZ 
states respectively. Pairs in which one fermion is outgoing and one incoming 
are associated with t-channcl enhancements. For each set we determine the 
maximum number of enhancements which can be simultaneously present and 
introduce the corresponding channel in the list of channels utilized for the 
process. Multi-mapping and adaptivity will take care of all related partially- 
resonant or non-resonant configurations. Each channel in Eq.(l) is integrated 
separately with VEGAS. 

In Fig. 9 we show pictorially the fiexible phase space mappings which are 
available in PHANTOM. We denote by Jl^ the usual angular variables in the 
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Fig. 9. Examples of phase space mappings 

center of mass of particles i and j and define Sj^...j^ = (pj^ + ... + Pi„)^ and 
tij = {Pi ~Pjy- Then, for instance, the mapping corresponding to the central 
figure of the first row will use as variables tn, t28, S3456, f^(34)(56)5 -534, ^^34, 
■556, ^56, ■534567- The remaining variables are the azimuthal angles 08 and 07 of 
particle 8 in the overall centre of mass and of particle 7 in the center of mass of 
34567, respectively. If we indicate with Mx the invariant mass at which the sx 
variable can resonate or the mass of the vector boson which enters into the tx 
channel exchange, we could have M17 = M28 = M34 = Msg = M^y, M345e = 
Mh, Mz which describes a WW fusion channel with a Higgs resonance; on the 
other hand we could also have M17 = ,M28 = M56 = My/, M34 = Mz which 
describes a '-)W ZW virtual scattering. It is clear that a large number of 
possibile channels can be described by this single parametrization. 



4 Running the program 

Modes of operation 

The program has two modes of operation which are selected by the input 
value ionesh. If ionesh=0 the program computes the cross section for one 
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single process, specified by iproc. If the input variable if lat is set to 1 the 
program produces one or more integration grids depending on the number of 
channels required for a good mapping of phase space for the selected process. 
If ionesh=l the program generates unweighted events for a set of processes, 
which is specified by the user, using the previously produced grids. All the 
integration grids for each selected process must be included for a meaningful 
generation. 

When ionesh=0 the program proceeds in two steps. The first one is called 
thermalization. It determines the relative weight of each channel in the mul- 
tichannel integration and it produces a first instance of PHANTOM space grids, 
one per channel. At least three thermalization iterations are performed. At 
the end of the third iteration all channel whose relative weight is smaller than 
10~^ are eliminated. If any channel is discarded, two additional thermalization 
iterations are performed, regardless of the fiag itmx.therm mentioned below. 
The grids produced in the thermalization stage are then used as a starting 
point for the second step which consists of one integration per channel. Each 
integration will typically consist of several iterations and at each iteration the 
PHANTOM space grid will be refined in an efi'ort to decrease the overall variance. 
A number of iterations between 3 and 5 is normally the best choice. If higher 
precision is requested it is usually more convenient to increase ncall rather 
than itmx. The user must be aware of the fact that if no point survives the 
cuts during an iteration, either during thermalization or at the integration 
stage, VEGAS will stop with an error. 

At pp and e'^e~ colliders, when ionesh=0 the process is computed exactly as 

specified by the user, assuming the first incoming particle to be moving in 
the +z direction and the second one in the —z direction. Even though at a 
pp collider this violates the symmetry between the two initial state protons, 
it can be useful for testing and for speciahzed studies. The full cross section 
can be easily obtained for unlike incoming partons by symmetrizing all dis- 
tributions with respect to the beams and multiplying by a factor of two. At 
a pp collider the default for unlike incoming partons is to sum over the two 
possible assignement of the two partons to the beams. The behaviour at pp 
colliders can, and typically should, be modified when ionesh=l by the fiag 
i.exchincoming. 



4-2 Input 

The syntax of the input is almost identical to the one required by the CERN li- 
brary routine FFREAD. Routines internal to PHANTOM are however used (iread, 
rread), so that real variables can (and must) be given in double precision. 
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All lines in the file of input must not exceed 80 characters, with the exception 
of filenames which can be 200 characters long. A * or C character at the 
beginning of a line identifies it as a comment line. Comment lines can be 
freely interspersed within the input file, with the only obvious exception that 
they must not interrupt a list of input values for a single array variable. The 
name of the variable to be read must be specified as the first word of a line 
(needs not to begin in column 1). Its valuc(s) must follow it. The list of values 
can span several lines. Variables which are not needed for the process under 
study will be ignored. They can be left in the file without harm. All variables 
actually read from the input file will be reproduced in the output. 

The input variable iproc specifies the desired process using the standard 
Monte Carlo particle numbering scheme: 



d 


u 


s 


c 


b 


e 








r 




9 


1 


2 


3 


4 


5 


11 


12 


13 


14 


15 


16 


21 



Antiparticle are coded with the opposite sign. The first two entries represent 
the initial state partons. Therefore the string 

3 -4 2 - 2 3 - 3 13 - 14 (3) 



corresponds to the reaction 

sc — > uussiiV^. (4) 



4.3 Cuts 



The program provides two means of specifying acceptance cuts. A basic pre- 
defined set can be imposed through the input file. They cover the usual re- 
quirement of minimum energy, transverse momentum and separation among 
the final state partons. All cuts are specified by an integer of the type i_f lag 
which specifies whether the corresponding cut is activated (i_f lag=l) or not 
(i_f lag=0) and by one or more values which define the extrema of the accepted 
region. The name of the flag in most cases is the name of the corresponding 
variable with i_ prepended. 

The flrst part of the input cuts is common to both running modes, and must 
be always kept unchanged when passing from ionesh=0 to ionesh=l. It con- 
stitutes in fact the setup under which phase-space grids are produced. In order 
to give the possibility of imposing other cuts at generation level, the input- 
file has also a cut section specific of the one-shot mode. The corresponding 
variables are the same as those in the common input section; they are just 
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renamed with a suffix - os - appended. These additional cuts, operating at 
generation level, are obviously effective only if more restrictive than the com- 
mon ones. Since no predefined set of cuts will be able to cover all possible user 
needs, it is also possible to include extra user-specified cuts via a routine called 
iuserfunc, an example of which is provided in the program package, inside 
the file cuts.f As for the predifined set of cuts, additional requirements can 
be imposed at generation level via the iuserf uncos routine. If user-specified 
cuts are imposed the package needs to be recompiled. 

4-4 Output 

The output of the program depends obviously on the mode of operation. The 
result of running the program in the integration phase (ionesh=0) is an output 

file in which are reported the input choices and the result of the various itera- 
tions during integration for the different channels. One can check from this file 
the accuracy reached by the integration and the reliability of the result. The 
informations needed for the successive unweighted event generations are given 
for channel xx of a given process by files under the name phavegasxx.dat. 
For every process a separate directory must be created. The output of the 
(ionesh=l) mode is a file in which the statistics for each channel contributing 
to the generation are reported, as well as the total integral of all channels 
contributing to the generations. It is advisable to check that this result cor- 
responds, within the statistical errors, to the sum of the cross sections of the 
various processes contributing to the final sample, taking into account possi- 
ble factors of two due to the exchange of incoming particles. If one has set 
in input iwrite_event 1. as it is normally the case, the unweighted events 
generated at parton level are written in a file named phamom . dat using the 
Les Houches Accord File Format [45] . The way the events are written depends 
on the fiag iwritejnothers. If this is set to 1 additional information is added 
to the event concerning possible "mothers" (resonances) which generate final 
particles or intermediate "on shell" particles. This information is used by the 
hadronization Monte Carlo's and in particular by PYTHIA [46] for the evolu- 
tion of the parton shower and the hadronization. Preliminary studies show 
that the use of this information has sizeable consequences, at least for top 
processes, on the final fully hadronized events. ^ It is obvious that in a full 
calculation, the identification of the "mothers" in an event is not straight- 
forward as it would be in a production times decay approximation. We have 
therefore chosen to somehow rely on the kinematic configuration of the event 
for assigning the "mothers" . We first of all determine which are the possible 
mothers considering the particular phase space channel and parameters used 



^ We thank Roberto Chierici for pointing out this feature to us and Torbjorn Sjos- 
trand and Fabio Maltoni for discussions on this point. 
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to produce the kinematical configuration. Wc tiicn look at the color configu- 
ration chosen via the amplitudes and test if the two informations agree. Only 
to those resonances which pass the check we assign the role of "mothers" . 



4.5 Scripts 



Determining the minimal set of reactions for which grids need to be computed 
can be tedious and error prone. Therefore we have included in the distribution 
two Perl scripts setupdir-all-LSF.pl and setupdir-all-LSF_ILC.pl which 
handle the task for hadron and e'^e~ colliders respectively. The script can be 
executed with, for instance: 

perl setupdir-all-LSF.pl [-options ...] 



where options include: 



-basedir 

-dirtreeroot 
-template 
-executable 
-inputstring 



-quarks 
-Top 



[./] 

[mhO] 

[template . stO] 
[phantom. exe] 
[""] 



-help 



directory which contains the exe file 

(full pathname) 
root of new tree (full pathname) 
input template file 
executable file 

must contain only leptons and gluons 
must be enclosed in double quotes : 

"e e_ mu vm_" 
only processes containing -inputstring 
will be generated 
number of quarks to be inserted 

(even integer > and <= 10) 
number of intermediate top/antitop quarks which can 
be reconstructed by the final state particles 
(possible values 0,1,2). If the option is of the form 
n+ any reaction with at least n intermediate tops are 
accepted. If the option is not set any number of tops 
is accepted (equivalent to 0+) . 
prints usage details 



[0] 



Long options can be abbreviated up to one letter, eg. -basedir X can be 
passed as -b X. In square brackets the default values are reported. For each 
reaction the scripts create a directory which contains the corresponding r . in 
input file which is generated using the given template which must contain all 
input parameters with the exception of iproc which is supplied by the script. 
A runfile called run is also created. In the root directory of the new tree 
the script creates a file called LSFf ile which can be used for submitting all 
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integration jobs to the CERN LSF batch system. For different batch systems 
the user should appropriately edit the Perl scripts. 



5 0{a%j^)-\-0{a'^j^jal) results in the semileptonic //z/^ channel 

In this section we present an example of some phenomenological results ob- 
tained with PHANTOM. These studies have been performed via the generation 
of high statistic samples of unweighted events at parton level for all processes 
with a pair of leptons iiu^ and different Higgs masses at the LHC. For a more 
complete discussion of other studies performed with PHANTOM we refer to [47] . 

We investigate how much the sensitivity to the VV scattering signal in the 
semileptonic /xi/^ channel is affected by the QCD irreducible background. 

In the absence of firm predictions in the strong scattering regime, trying to 
gauge the possibilities of discovering signals of new physics at the LHC requires 
the somewhat arbitrary definition of a model of V^Vl scattering beyond the 
boundaries of the SM. Some of these models predict the formation of spectac- 
ular resonances which will be easily detected. For some other set of parameters 
in the models only rather small effects are expected [48,49]. 

Two groups [50,51] have recently tried to parametrize the low energy behaviour 
of large classes of composite Higgs models. They showed that the most char- 
acteristic signature is a reduced coupling between the Higgs state and the SM 
vector bosons, leading to only a partial cancellation of the growth of Yl^l 
scattering amplitudes which would eventually result in violations of unitarity. 

The predictions of SM in the presence of a very heavy Higgs provide an up- 
per bound on the observability of such effects. The linear rise of the cross 
section with the invariant mass squared [52] in the hard VV scattering will 
be swamped by the decrease of the parton luminosities at large momentum 
fractions and, as a consequence, will be particularly challenging to detect. 

In our study, we consider the full set of parton-level processes involved at 

qq qqqqi^v^, gg qqqqi^v^ 
gq gqqqi^P„ qq ^ ggqqiJ^v^ 

together with a selection procedure as close as possible to the actual experi- 
mental practice, without resorting to any flavour information other than the 
one which a typical 6-tagging algorithm is able to provide. A more complete 
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analysis including the AjW background at 0{a1j^ag) is left for a forthcoming 
paper. 

The Pz of the neutrino is approximately reconstructed requiring the invari- 
ant mass of the two leptons to be equal to the W boson nominal mass. Its 
transverse momentum is assumed to be equal to the missing px- 

All the results presented in this section have been obtained using the CTEQ5L 
PDF[53] set. The QCD couphng constant has been evaluated at the scale 

Q' = M^^ + pT{topr, (5) 



where Pritop) is the transverse momentum of the reconstructed top, for all 
processes in which a t or i can be produced. For all other processes the scale 
has been evaluated as 

Q'-M^+lj:pi,, (6) 



where pn denotes the transwrso momontum of the ?-th final state particle. 



Mh = 200 GeV 




(^EW+QCD 


all events 


0.89 pb 


80.8 pb 


top events 


0.52 pb 


71.6 pb 


ratio top/all 


0.58 


0.89 



Table 1 

Contribution of tt/ single t to the total cross section with standard acceptance cuts 
only (see the left part of Tab. 2). Comparison between results at 0{a%j^f) (EW) 
and C'(a|^^)+C'(a|^a|) (EW+QCD). Interferences between the two perturbative 
orders are neglected. 

It should be clear from the results shown in Tab.l that suppressing the top 
background is the primary goal to achieve. In this analysis we assume the 
possibility to tag 6-jets in the region < 1.5 with 0.8 efficiency, which allows 
to discard part of the events containing b quarks in the final state. We impose 
additional cuts against the top on the invariant mass of triplets of type {jjj} 
and {jfiu}, where j denotes any final-state quark or gluon. In order to isolate 
two vector boson production, kinematical cuts are applied on the invariant 
mass of the two most central jets, which are associated in our analysis to a 
W or Z decaying hadronically. The VV fusion signature is further isolated by 
requiring a minimum Ar) separation between the two forward/backward jets. 

At this stage, however, any attempt to appreciate differences between Higgs 
and no- Higgs scenarios at high invariant masses would still be vain. This is 
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Acceptance cuts 

Prii^J) > 10 GeV 



E{£^,j) > 20 GeV 



<6.5 



M{jj) > 20 GeV 



Selection cuts 

6-tagging for |r/| < 1.5 (80% efficiency) 



\M{jjj;j£^u, 



Mtop\ > 15 GeV 



70 GeV < M(jcjc) < 100 GeV 
M{jfjb) < 70 GcV ; M{jfjb) > 100 GeV 



M{jj) > 60 GeV 



PT(^^Z^rec) > 100 GeV 



M{jf/i,e^urec) > 250 GeV 

Table 2 

List of kinematical cuts applied in all results on the /iv^ channel, j denotes any 
final-state quark or gluon, while £^ is the charged lepton. The subfixes c,f,b mean 
central, forward, backward respectively, i^rec is the neutrino reconstructed following 
the prescription described in the text. 





no Higgs 


Mh = 200 GcV 




a 


£=100 fb"^ 


a 


£=100 fb"^ 


ratio 


all events 


12.46 fb 


1246 ± 35 


13.57 fb 


1357 ± 37 


0.918 


Me„t = 0.8 TeV 


3.19 fb 


319 ± 18 


1.45 fb 


145 ± 12 


2.200 


M,ut = 1-2 TeV 


1.28 fb 


128 ± 11 


0.41 fb 


41 ± 6 


3.122 


Mcut = 1-6 TeV 


0.60 fb 


60 ± 8 


0.14 fb 


14 ± 4 


4.286 


le 3 



Integrated 0{a%j^) cross section for M{jcjcli^) > Mcut and number of expected 
events after one year at high luminosity having applied the cuts listed in Tab. 2. 





no Higgs 


Mh = 200 GeV 




a 


£=100 fb-^ 


a 


£=100 fb-^ 


ratio 


all events 


40.70 fb 


4070 ± 64 


40.73 fb 


4073 ± 64 


0.999 


Mcut = 0.8 TcV 


7.61 fb 


761 ± 28 


5.14 fb 


514 ± 23 


1.481 


Mcut = 1-2 TeV 


2.53 fb 


253 lb 16 


1.73 fb 


173 it 13 


1.462 


Mcut = 1.6 TeV 


1.00 fb 


100 ± 10 


0.55 fb 


55 ± 7 


1.818 



Table 4 

Integrated 0{a%j^f) + 0{a%j^a1) cross section for M{jcjcl^) > Meat and number of 
expected events after one year at high luminosity having applied the cuts listed in 
Tab. 2. Interferences between the two perturbative orders are neglected. 
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essentially due to the fact that the contribution of the QCD diagrams is not 
substantially affected by the above-mentioned selection criteria. Investigating 
the differences between the kinematics of VV scattering and VV + 2 jets, we 
have identified additional cuts that serve our purpose. As the background 
dominates in the phase space regions characterized by vector boson emitted 
by the tag jets, a viable method of taming VV + 2 jets consists in applying 
cuts on the pt and rj of the W reconstructed from leptons as well as on the 
invariant mass of the W plus any of the two tag jets. All details about the 
selection cuts applied are reported in Tab. 2. 



Mh = 200 GeV I No Higgs 




M||„ (GeV) (GeV) 



Fig. 10. Invariant mass distribution of the two leptons and the two most central jets 
in the Standard Model with a light Higgs (on the left) and in the no-Higgs scenario 
(on the right). The cuts applied are listed in Tab.2. 0{a%,,) (EW) and ©(ol^al) 
(QCD) contributions to the differential cross section have been isolated and are 
shown separately. The QCD contributions are further split into top background (in 
blue) and VV + 2 jets (in green). 

Fig. 10 illustrates the final results of our analysis, showing that the top back- 
ground is basically under control. VV + 2 jets still provides a non-negligible 
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contribution over the whole invariant mass spectrum, nevertheless differences 
between light-Higgs and no-Higgs can be appreciated. 

In Tab. 3,4 we show the integrated cross section at high energies as a function 
of the minimum invariant mass, comparing results for the pure EW and the 
EW+QCD cases. Despite reducing the ratio between no-Higgs and light-Higgs 
cross sections, the inclusion of QCD background does not seem to compromise 
the possibility of detecting signals of EWSB at the LHC from an excess of scat- 
tering events. We find that about 500 events are expected above 800 GeV after 
one year of high luminosity running {C = 100 fb~^) in case of a Higgs boson 
with mass 200 GeV. The no-Higgs model predicts about 250 more events in 
accordance with the enhancement of the VV differential cross section at high 
energies. These numbers refer to the muon channel only, and can obviously 
be improved by summing up the muon and electron channels. It should nev- 
ertheless be noticed that imposing a minimum AR separation among colored 
particles could degrade these preliminary results and requires further investi- 
gations. 



6 Conclusions 

We have described in detail the features of the Monte Carlo PHANTOM which 
is the first dedicated event generator for six fermion physics which can be 
used at high energy pp, pp and e'^e~ colliders. It has already been used for 
some phenomenological studies on boson boson scattering physics at the LHC 
and a short example of its possibilities in this field has been presented. It 
will be used for further detailed analyses at parton level, which should also 
contribute to assess the complementary role of hadronic and e~^e~ colliders in 
this field. It is however important to stress that after partonic studies, more 
realistic ones with full event reconstruction and detector simulation are needed 
to fully exploit the physics possibihties. PHANTOM has been developed also as 
a tool for these more complete analyses. 

The program can be downloaded from http: / /www.ph. unito.it/~ballestr/phantom 
where all new versions will be made available. 
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Appendices 



A Routines 

Steering routines. The first one is the main program, the second drives the 
integration. The third, in generation mode, reads all grid files and steers the 
unweighted sampling. 

phantom. f integ.f oneshot.f 

Routine which parses the input file, 
readinput . f 

As the name indicates coupling. f computes the EW couplings in the G^- 
scheme, set masses and computes decay widths. 

coupling. f 

Modified VEGAS package. 

phavegas . f 
Integration function for phavegas. f 

fxn.f 

These routines select the set of channels to be used for each process and 
initialize the corresponding parameters. 

proc.f procini.f procextraini .f 

Phase space routines. The XYZ j ac . f routines, not shown here, compute the 
jacobian for the corresponding channel from the momenta. 

phspl_l_31_multi_c . f phspl_l_4_multi5_c . f 

phspl_2_3_multi5_c . f phspl_5tol_4to31_multi_c . f 

phsp2_4_multi5 . f phsp2_4to31_multi5 . f 
phsp3_3_multi5 . f 

Beamstrahlung[54] and Initial State Radiation routines for e+e~ colliders. 

circe.f isr.f 

Implementing the cuts 
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cuts . f 



These routines compute the complete amplitude for each process using the 
master amplitudes and the color evaluation routines mentioned below. 



Color evaluation routines 

colevalew.f coleval.f 
Storing the event information according to the latest Les Houches Accord 

LHAFilelnit.f storeLH.f 

Utility routines 

phread.f bernoul.f util.f ccfcsym.f 
extrema.f isign.f perm2g.f perm.f 



B Sample input file 

We report here a sample imput file to integrate the process bb ^bb fi'^fi~ p^. 
The hues starting with * and the text foUowing a ! are comments which can be 
left without problem in the input file. One can use this same file to generate 
an unweighted sample with the same cuts and options, just changing ionesh 
in ionesh 1 and fixing at the end the exact number of files nf iles and 
their names for the processes to be generated. The line containing the indica- 
tion of the process (iproc) will in this case be ignored. It is evident from the 
comments that some parts of the input arc ignored when ionesh=0 and some 
others when ionesh=l. We also recall that our convention is that whenever 
there is a yes/no option 1 corresponds to yes and to no. The EVENTUAL MORE 
RESTRICTIVE CUTS are not all reported for brevity. 

iproc -5 -5 -5 -5 13 -13 14 -14 

idum -123456789 ! random number seed must be a large negative number 



ampem . f ampSf qcd . f amp2g . f 



Master amplitudes: 4W, 4Z, 2W2Z, 2g3Z, ggW2Z. 



fourw.f fourz.f 
fourwqcd.f fourzqcd.f 
ggSz.f ggzww.f 



twoztwow . f 
twoztwowqcd. f 
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PDFname /home/phantom/phantom/lhapdf -5 . 2 . 2/PDFsets/cteq51 . LHgrid 



* i_PDFscale selects the PDF scale: 

! =1 for all processes, based on pT's of ALL OUTGOING PARTICLES 
! =2 process by process, based on pT of the (RECONSTRUCTED) TOP 
! if possible, otherwise as done in option 1 

i_PDFscale 2 



* i_coll determines the type of accelerator: 

* 1 => p-p 2 => p-pbar 3 => e+e- 
i_coll 1 

i_isr ! yes/no ISR for e+e- collider only 

i_beamstrahlung ! yes/no beamstrahlung for e+e- collider only 

* perturbativeorder = 1 alpha_em~6 with dedicated amp 

* 2 alpha_s''2alpha_em~4 3 alpha_em~6 + alpha_s"2alpha_em"4 
perturbativeorder 3 

ionesh ! 0= one process run, 1= one shot generation 
ecoll 14000. dO ! collider energy 



rmh -500. dO ! Higgs mass (GeV) . If negative Higgs diagrams are not 

! computed: equivalent to rmh=inf inity . 

i_ccfam 1 ! yes/no family+CC conjugate 

* CUTS 



i_e_min_lep 1 
e_min_lep 20 . dO 

i_pt_min_lep 1 
pt_min_lep 10. dO 

i_eta_max_onelep 

eta_max_onelep 3 . dO 

i_eta_max_lep 1 
eta_max_lep 3 . dO 

i_ptmiss_min 
ptmiss_min 50. dO 



! yes/no lepton energy lower cuts (GeV) 
! yes/no lepton pt lower cuts (GeV) 



yes/no at least one lepton in absval of 
eta_max_onelep . If no leptons are present 
in the final state cut is ignored 



! yes/no lepton rapidity upper cuts 
! yes/no missing pt lower cuts 
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i_e_niin_j 
e_miii_j 

i_pt_min_j 
pt_min_j 



1 

20. dO 

1 

10. dO 



! yes/no jet energy lower cuts (GeV) 
! yes/no jet pt lower cuts (GeV) 



i_eta_max_j 1 ! yes/no jet rapidity upper cuts 

eta_max_ j 6 . 5d0 



i_eta_jf _jb_jc 
eta_def _jf _min l.dO 
eta_def _jb_max -l.dO 
eta_def _jc_inax 3.d0 



yes/no the following 3 cuts 
min rapidity for a jet to be called forward 
max rapidity for a jet to be called backward 
max rapidity for a jet to be called central 
(absval) 



i_pt_min_jcjc 
pt_min_jcjc 50. dO 

i_rm_min_jj 1 
rm_min_ j j 20 . dO 



! pt lower cut on two central jets (GeV) 
! yes/no minimum invariant mass between jets (GeV) 



i_rm_min_ll 1 ! yes/no minimum invariant mass between charged lept. 

rm_min_ll 20. dO 

i_rm_min_jlep ! yes/no minimum invariant mass between jets and lepton 
rm_min_jlep 30. dO 

i_rm_min_jcjc ! yes/no minimvim invariant mass between central jets 

rm_min_jcjc 20. dO 

i_rm_max_jcjc ! yes/no maximum invariant mass between central jets 
rm_max_jcjc 14000. dO 

i_rm_min_jf jb ! yes/no minimum invariant mass between forward 

land backward jets 

rm_min_jfjb 100. dO 

i_eta_min_jf jb ! yes/no minimum rapidity difference between forward 
! and backward jet 

eta_min_jfjb 3.d0 

i_d_ar_jj ! yes/no minimum delta_R separation between jets 
d_ar_jj 0.7d0 

i_d_ar_jlep ! yes/no minimum delta_R separation between jets and lepton 
d_ar_jlep 0.7d0 
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i_thetamin_j j ! yes/no mininnam angle separation between jets (degrees) 
thetamin_jj 15. dO 

i_thetaiiiin_jlep ! yes/no minimum angle separation between jets 

! and lepton (degrees) 
thetamin_jlep 15. dO 

i_thetaiiiin_leplep ! yes/no minimum angle separation between charged 

! leptons (cosine) 
thetamin_leplep 15. dO 

i_usercuts ! yes/no (1/0) additional user-defined cuts 



****** IF (lONESH.EQ.O) THEN 



acc therm 



O.OldO 



! thermalization accuracy 



ncall_therm 300000 3000000 ! thermalization calls per iteration 

! for the first 3 and for the remaining iterations 
itmx therm 5 ! thermalization iterations 



acc 
ncall 
itmx 
iflat 



0.005d0 ! integration accuracy 

5000000 ! integration calls per iteration 

5 ! integration iterations 

1 ! yes/no flat event generation 



****** ELSEIF (lONESH.EQ.l) THEN 



scalemax 1 . IdO 



! scale factor for the maximum 



nunwevts 20000 ! number of imweighted events to be produced 

iwrite_event 1 ! yes/no momenta of flat events written in .dat files 

iwrite_mothers 1 ! yes/no information about intermediate 

! particles (mothers) in .dat files 

i_exchincoming 1 



EVENTUAL MORE RESTRICTIVE CUTS: 

SAME LIST AS BEFORE WITH "os" POSTPENDED 
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iextracuts ! yes/no more restrictive cuts at generation level 

i_e_iiiin_lepos ! lepton energy lower cuts (GeV) 

e_min_lepos 30. dO 

i_usercutsos 



* nfiles= number of files from which take the input for 

* one shot generation 

nfiles 3 ! number of input files for oneshot generation 

/home/phantom/phantom/processl/phavegasOl . dat 
/home/phantom/phantom/processl/phavegas02 . dat 
/home/phantom/phantom/process2/phavegas01 . dat 
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